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Abstract 

•  For  a  fev  particular  vavegulde  problems^  standard  variational 
expressions  have  pre\T.ously  been  shown  to  be  upper  or  lower  bounds  on 
the  quantities  of  interest.   However,  bounds  have  not  previously  been 
obtained  for  any  truly  three  dimensional  problem,  that  is,  where  the 
fields  cannot  be  derived  from  a  single  scalar  potential.  An  example  is 
a  three  dimensional  obstacle  which  contacts  only  one  waveguide  surface. 
As  one  consequence,  no  straightforward  procedure  exists  for  improving 
the  approximations.   Recently,  Kato  devised  a  rather  general  method  for 
bounding  the  cotangent  of  the  phase  shift  for  a  given  angular  momentum 
in  a  quantum  mechanical  central  potential  scattering  problem.   This 
method  should  be  applicable  whenever  a  system  can  be  analyzed  in  terms 
of  uncoupled  standing  waves  each  characterizable  by  one  real  phase  shift. 
The  method  is  here  adapted  to  waveguides,  including  truly  three  dimen- 
sional problems.   The  obstacle  must  be  symmetric  about  a  plane  perpen- 
dicular to  the  waveguide  axis,  with  certain  exceptions  only  one  mode 
should  propagate,  and  the  system  should  be  lossless.  AnaJLysis  is  then 
possible  in  terms  of  i]  •  and  r|  ,  the  real  uncoupled  phase  shifts  associated 
with  the  even  and  odd  standing  waves,  respectively.   The  bounds  obtained 
on  cott]  and  cotn  determine  bounds  on  the  equivalent  network  elements. 
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1.   Introduction 

The  introduction  of  variational  techniques  has  greatly  facilitated 
the  solution  of  scattering  problems  in  various  fields.   They  were  origin- 
ally applied  by  Schwinger  to  scattering  in  an  analysis  of  waveguide 

If  2  3 

problems   '   and  quantum  mechanical  problems   and,  independently. 

X  k  5 

by  Hulthen   in  an  analysis  of  quantum  mechanical  problems.    These 

techniques  have  subsequently  proved  to  be  useful  in  a  wide  variety  of 
fields  including  optics,  acoustics  and  water  waves.   Though  they  are  clearly 
superior  to  previously-used  techniques,  variational  techniques   as  applied 
to  scattering  theory  suffer  a  major  defect  as  compared  to  variational  tech- 
niques as  applied  to  the  determination  of  point  eigenvalues  such  as  the 
resonant  frequencies  of  vibrating  systems  or  the  bound  state  energy  eigen- 
values of  quantum  mechanical  systems.   This  defect  is  that  while  the  expres- 
sions are  stationary,  they  are  in  general  neither  an  upper  nor  a  lower  bound. 
Thus,  given  two  different  trial  functions,  one  cannot  generally  know  which  of 
the  two  variational  results  is  better;   in  particular,  the  incorporation  of 
additional  parameters  in  a  given  trial  function  does  not  guarantee  a  better 
result.   Consequently,  it  would  be  desirable  to  be  able  to  deduce  bounds  on 
the  quantities  of  interest.   Now,  in  fact  in  some  particular  waveguide  pre- 

blems .  the  variational  expressions  themselves  have  been  shown  to  be  upper  or 

1.2 
lower  bounds.      These  problems  are  not  sufficiently  general,  however. 

to  include  all  cases  of  practical  interest.   For  example,  no  technique  seems 
to  have  been  devised  which  is  capable  of  giving  bounds  on  the  relevant  quan- 
tities for  a  truly  three-dimensional  case,  that  is,  where  the  fields  cannot 
be  derived  from  a  single  scalar  potential. 
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Nov  a  rather  general  approach  has  been  devised  by  T.  Kato      ' 
for  obtaining  boimds  on  cot  (r|  -  O),  where  r)  is  the  phase  shift  for  a  given 
angular  momentum  in  a  quantum  mechanical  central  potential  scattering  problem 
and  9  is  a  parameter  chosen  for  convenience.   This  method  should  be  applicable 
to  any  problem  which  can  be  analysed  in  terms  of  uncoupled  standing  waves  each 
of  which  can  be  characterized  by  one  real  phase  shift.   It  is  the  purpose  of 
this  paper  to  adapt  the  Kato  fonnalism  to  scattering  by  obstacles  in  wave- 
guides, including  the  truly  three  dimensional  problem.   The  problems  to  be 
considered,  for  reasons  which  will  be  discussed  later,  are  subject  to  the 
following  restrictions: 

(1)  The  waveguide  is  uniform  in  the  asymptotic  region. 

(2)  Only  the  dominant  mode  propagates,  with  exceptions  which  will  be 
mentioned  later. 

(5)    Tlie  obstacle  is  symmetrical  about  a  plane  perpendicular  to  the 
axis  of  the  waveguide. 

{h)        The  waveguide  and  obstacle  are  lossless. 
In  the  companion  paper.     the  Kato  formalism  is  applied  to  some  specific 
problems  involving  dielectric  obstacles  and  it  is  found  that  with  a  simple 
trial  function  and  a  moderate  amount  of  work  one  can  obtain  upper  and  lower 
bounds  on  cot  (t|  -  O)  and  cot  {r\     -  0)  where  t),  and  r]  are  the  real  uncoupled 
phase  shifts  associated  with  the  even  and  odd  standing  waves,  respectively. 
These  bounds  differ  from  one  another  by  only  a  few  per  cent,  even  for  obstacles 
whose  dimensions  are  comparable  with  the  transverse  dimensions  of  the  waveguide. 
In  addition,  for  the  obstacles  considered,  bounds  can  be  obtained  directly  on 
the  phase  shifts  themselves  and  on  the  elements  of  the  equivalent  'T'  network. 
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2.   Applicability  of  the  Kato  method 

Before  proceeding  to  the  waveguide  case^  it  will  he  helpful  to  outline 
Kato's  approach.   A  standard  technique  in  the  quantiim  theory  of  scattering 
of  a  particle  by  a  central  potential,  V(r)  =  V(r) ^  is  to  separate  the  wave 
into  partial  waves,  i.e.,  the  eigenfunctions  of  the  angular  momentum  operator. 
For  a  specified  angular  momentum  and  energy,  the  effect  of  the  potential  on 
each  incident  partial  wave  is  completely  characterized  by  one  real  number, 
the  phase  shift.   For  the  case  of  zero  orbital  angular  momentum,  which  is  the 
closest  analogue  to  the  waveguide  case  we  are  ultimately  interested  in,  the 
phase  shift  t|  is  determined  by  the  equations 


j:\  = 


^  +  k^  +  W(r) 


dr^ 


^^Q  =  0   '   0  <  r  <  oo        (la) 


0 


Uq(0)  =  0  (lb) 


u  (r)  -^  cos(kr-i-©)  +  cot(r|-0)  sin(kr+0)   ,  r  -^  oo   (ic) 


2 
In  the  above  expressions,  k  and  W(r)  are  proportional  to  the  total  energy  and 

the  negative  of  the  potential  energy,  respectively,  and  the  normalization  con- 
stant 0  satisfies  0  <  ©  <  n  but  is  otherwise  arbitrary.   The  introduction  of 
0  provides  an  additional  element  of  freedom. 

Let  u  be  a  trial  wave  function  which  satisfies  the  boundary  condition 
(lb)  and  has  the  asymptotic  form  (ic)  with  a  trial  phase  shift  t]  ,  but  which 
does  not  satisfy  the  wave  equation  (la).  (Trial  quantities  will  be  distinguished 


h   - 


from  their  exact  counterparts  by  a  subscript  t.)   From  Green's  theorem. 


I 


/  /         t-     du^,       du„-i 


9t       

^9  dr  "  ^9t  dr 


o 


(2) 


and  the  asymptotic  forms  of  u  and  u  ,  it  follows  that 


k  cot(Ti-9)  =  k  cot(Ti^-9)  -  j  u^^aC  u^^dr  +  j  Wg<^  Wq  dr  ,      (3; 


where 


vjr)   =  u^Jr)  -  u  rr)   .  (1^.) 


'9^-'      Gf '    0 


Equation  (5)  is  an  identity.   If  u   is  a  good  approximation  to  u^,  then  w^ 


is  a  small  term  of  first  order,  and  the  error  tenn   (  v^<^  w  dr  is  of  second 


j  w  <3^  w 

order.   If  the  error  term  is  neglected  in  equation  (2),  the  remaining  two 
terms  on  the  right  hand  side  can  be  evaluated  explicitly  and  constitute  a 
variational  expression  for  k  cot(ri-9).   (For  certain  values  of  9,  this 
expression  reduces  to  previously  known  variational  principles.) 

Kato's  major  contribution  was  the  determination  of  upper  and  lower  bounis 
on  the  error  term  and,  hence,  on  k  cot(T|-9).     Kato  derives  the  inequality 


%'^  j  p'^U'^e^f^^  £  J  Vq^'V^  S  Pg"^/  P'^ic^'n^^fdv 


(5) 
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where  p  is  some  non-negative  weight  factor  to  be  chosen  for  convenience; 
a  is  the  smallest  positive  eigenvalue  and  p  the  smallest  (in  absolute 
value)  negative  eigenvalue  of  the  associated  eigenvalue  problem, 

o<^  0  (r)  +  [1  p(r)  0  (r)  =  0   ,   0  <  r  <  oo  .  (6) 


The  eigenfunctions  p  and  eigenvalues  p.  are  determined  by  the  boundary 
conditions  that  0  vanish  at  the  origin  and  have  the  asymptotic  phase  shift, 
5(n^)j  given  by 


6(n  )  =  0  +  nrt   ,   n  =  0,  +  1, 


n 


(T) 


Even  quite  crude  lower  bounds  on  a^  and  p  in  equation  (5)  will  provide  close 
bounds  on  the  error  term  in  equation  (3)  and,  hence,  on  the  phase  shift,  pro- 


-,    j  p"  {Xn 


vided  the  error  integral,   |  p  ""(c^Uq,  )   dr,  which  vanishes  for  the  exact  wave 


function,  can  be  made  sufficiently  small.   Comparison  potentials  for  which 

equation  (6)  can  be  solved  exactly  are  employed  in  the  estimation  of  a_  and 

P  ,  together  with  a  theorem  which  states  that  the  phase  shift  increases 

6,7,8 
monotonically  with  increasing  'potential',  W(r)   (decreasing  potential  energy) 


Let  us  now  consider  a  one  dimensional  quantum  mechanical  scattering  pro- 
blem, which  is  a  still  closer  analogue  to  the  waveguide  problem.   This  differs 
from  the  previous  case  in  that  the  range  of  the  independent  variable  extends 
from  -  00  to  +  00  .   The  difference  is  fundamental,  however,  since  there  are 
now  two  channels,  i.e.,  two  sets  of  incoming  and  outgoing  waves.   In  a  scatter- 
ing problem  which  contains  n  channels,  the  remote  effects  of  the  scattering 
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process  are  characterized  by  an  n  x  n  scattering  matrix  relating  the  ampli- 
tudes of  the  n  outgoing  waves  to  those  of  the  n  incoming  vaves.   Since  the 

2  2 

n  elements  of  the  scattering  matrix  are  complex,  there  are  2n  quantities  in 

all.   However,  not  all  of  these  quantities  are  independent.   It  follows  from 
conservation  of  probability  (energy  in  the  electromagnetic  case)  and  time 
reversibility  that  the  scattering  matrix  is  unitary  and  symmetric,  which 
reduces  the  number  of  independent  quantities  to  1  n(n  +  l) .    "    In  the 
case  of  scattering  by  a  central  potential,  there  is  a  channel  corresponding 
to  each  value  of  the  angular  momentum,  but  since  the  potential  does  not  couple 
the  waves  of  different  angular  momentum  there  is  a  separate  scattering  problem 
for  each  channel  characterized  by  a  single  number,  the  phase  shift.   In  the 
one  dimensional  problem  there  are  two  channels;  hence  three  independent  quan- 
tities are  required  to  characterize  the  scattering.   For  this  case,  -the  Kato 
method  does  not  appear  to  be  applicable.   However,  if  we  restrict  ourselves  to 
an  even  'potential',  W(x)  =  W(-x),  the  number  of  independent  parameters  is 
reduced  to  two.   Further,  the  two  independent  solutions  of  the  wave  equation 
can  be  taken  as  standing  waves  which  are  even  ajid  odd  in  x,  each  characterized 
by  one  real  number,  the  phase  shifts  ti  and  r\   ,  respectively.   The  even  and 
odd  solutions  can  be  regarded  as  solutions  to  two  completely  independent 
scattering  problems,  distinguished  by  the  boundary  conditions.   The  odd  solu- 
tion is  identical  with  the  solution  of  the  radial  wave  equation  for  zero 
angular  momentum  (an  odd  function  vanishes  at  the  origin)  except  that  the 
integrals  are  multiplied  by  j  because  the  range  of  integration  has  been 
doubled.   The  Kato  method  can  also  be  applied  to  the  even  solution  with  minor 
modifications.   The  details  will  not  be  presented  here  because  of  the  close 
correspondence  to  the  waveguide  case  which  will  now  be  treated. 


In  the  waveguide  problem  two  additional  formal  complications  arise. 
Firstly,  the  problem  is  three  dimensional  so  that  the  line  integrals  which 
appear  in  the  variational  principle  ai«  replaced  by  volume  integrals.   Secondly, 
the  wave  functions  are  vectors  rather  than  scalars.   It  will  be  seen  that  these 
two  complications  do  not  cause  any  difficulty  in  principle;  in  particular,  the 
assumed  uniformity  of  the  waveguide  in  the  asymptotic  region  makes  the  problem 
effectively  one  dimensional  in  this  region.   However,  there  are  two  channels 
corresponding  to  each  propagating  mode.   If  the  Kato  method  is  to  be  applicable, 
then  in  line  with  the  previous  discussion  the  waveguide  and  obstacle  must  be 
such  that  it  is  possible  to  reduce  the  scattering  problem  to  completely  inde- 
pendent problems  each  characterized  by  one  real  number.   The  number  of  channels 
can  be  reduced  to  two  by  requiring  that  the  waveguide  have  only  one  propagating 
mode,  which  is  not  a  serious  restriction  in  most  applications.   Then,  as  in  the 
one  dimensional  case,  the  number  of  independent  parameters  can  be  reduced  from 
three  to  two  and  the  scattering  problem  to  completely  independent  problems 
for  even  and  odd  solutions  if  the  obstacle  is  symmetrical  about  a  plane  perpen- 
dicular to  the  axis  of  the  waveguide. 

More  than  one  mode  can  be  allowed  to  propagate  if  the  obstacle  is  such 
that  it  does  not  couple  these  modes.   As  an  example,  consider  a  metallic  or 
uniform  dielectric  obstacle  in  a  rectangular  waveguide  in  the  form  of  a  right 
circular  cylinder  parallel  to  the  direction  of  polarization  of  the  dominant 
mode,  centered  in  the  guide  and  extending  between  the  conducting  boundaries. 
It  is  possible  to  select  the  dimensions  of  the  waveguide  in  such  a  way  that 

in  some  range  of  frequencies  there  are  only  three  propagating  modes,  usually 

12 
designated  TE   ;.  '^'^on  ^^'^  '^'^m  '  These  three  modes  are  not  coupled  to  one 


another  by  such  an  obstacle,  nor  are  the  even  and  odd  standing  wave  solutions 
corresponding  to  each  mode  coupled.  Thus  there  is  a  total  of  six  channels 
such  that  the  scattering  in  each  channel  is  a  separate  problem  characterized 
by  a  single  phase  shift.  A  simpler  example  is  that  of  a  dielectric  obstacle 
whose  permittivity  is  an  even  function  of  the  coordinate  parallel  to  the  axis  of 

the  waveguide  only.   The  problem  is  then  one  dimensional  and  we  can  consider 
an  arbitrary  number  of  propagating  modes,  for  the  various  channels  are  uncoupled 
ajid  each  is  characterized  by  a  single  phase  shift.   The  problem  of  multiple 
propagating  modes  will  not  be  considered  further.   It  will  henceforth  always 
be  understood  that  just  the  dominant  mode  propagates. 

In  some  cases  the  symmetry  of  i:he  obstacle  is  such  that  the  scattering 
can  be  described  in  terms  of  a  scalar  wave  function  rather  than  a  vector  wave 
function,  which  may  result  in  considerable  simplification.   For  example^  in  a 
rectangular  waveguide  the  dominant  mode  can  be  described  variously  as  a  TE  wave 
propagating  parallel  to  the  axis  of  the  waveguide,  as  a  TM  wave  propagating 
parallel  to  the  axis  of  polarization,  or  as  a  TE  wave  propagating  in  a  direction 
perpendicular  to  both  of  these  axes.  Any  finite  obstacle  destroys  uniformity 
in  the  direction  of  the  waveguide,  but  it  may  possess  sufficient  symmefcrv'  to 
preserve  uniformity  in  one  of  the  perpendicular  directions,  in  which  case  the 
TE  or  TM  description  still  applies.   Now  IE   and  TM  waves  can  be  derived  from 
Hertzian  electric  and  magnetic  vector  potentials,  respectively,  which  are 
parallel  to  the  propagation  direction  of  the  TE  and  TM  waves  resx:iectlvely. 
Thus  the  fields  are  characterized  by  a  single  scalar  function,  the  component 
along  the  appropriate  axis  of  propagation  of  the  Hertzian  vector  potential. 


If  the  symmetry  of  the  obstacle  permits  description  of  the  fields  in  terms  of 
a  TM  wave  propagating  parallel  to  the  axis  of  polarization,  the  electric  field 
is  proportional  to  the  vector  potenti^il,  and  the  vector  equations  derived  below 
go  over  to  scalar  equations  directly.   For  the  symmetry  which  permits  descrip- 
tion in  terms  of  a  TE  wave,  there  is  no  such  simple  correspondence  and  the 
scalar  formulation  of  the  problem  requires  a  separate  derivation  starting 
directly  from  the  Hertzian  vector  potential  formulation.   Only  the  vector 
formulation  will  be  presented  here,  though  in  the  TE  case  it  might  be  simpler 
to  use  the  scalar  fonnulation. 

5.   Variational  principles  for  waveguides 

In  this  section  we  shall  derive  the  variational  princijrle  for  a  waveguide 
for  both  metallic  and  dielectric  obstacles,  with  the  four  restrictions  noted 
in  the  Introduction.   In  line  with  the  previous  discussion,  it  then  follows 
that  the  general  scattering  problem  reduces  to  two  independent  scattering 
problems,  each  of  which  can  be  characterized  by  just  one  real  phase  shift. 

Assuming  a  time  dependence  exp(-icDt),  the  electric  field  intensity, E(r) , 
satisfies  the  equations 

«^^  =  -  S.X  Y.X  ^  +   {^  1^   +  W)E  =  0  (8a) 

^  •  E.  =  0  (8b) 

where  c  is  the  velocity  of  light,  e  is  the  relative  permittivity,  and  where 

2        2 
the  'potential'  W  is  defined  by  W  =  cd  (e-  l)/c  ;  the  quantum  mechanical  analogue 

of  W  is  proportional  to  the  negative  of  the  potential.   In  addition,  E^  is 

required  to  be  normal  to  all  conducting  surfaces  including  any  metallic  obstacles. 
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Let  the  z-coordinate  be  parallel  to  the  axis  of  the  vaveguide.   It  follovs 
from  restriction  (3)  that  W(z)  =  W(-z),,  and  the  general  solution  can  be 
expressed  as  a  linear  comb-ination  of  standing  wave  solutions  which  are  even 
and  odd  functions  of  z.   In  what  follows,  subscripts  e  and  o  will  be  used 
to  denote  quantities  associated  with  the  even  and  odd  functions  respectively, 
but  will  be  omitted  wherever  the  formalism  is  the  same  for  both  functions. 
The  solutions  are  required  to  have  the  asymptotic  forms  for  z  -^  +  co 

E  (r)  =  e(x,y)  |-sin(kz+Q)  +  cot(n^-9)  cos(kz+0) |  (9a) 

E  (r)  =  e(x,y)  |  cos(kz+0)  +  cot(ri^-0)  sin(kz+Q)  |  (9b) 

where  e(x,y)  is  the  form  function  for  the  dominant  mode  and  where  0  satisfies 

0  <  9  <  jt.   The  utility  of  the  parameter  0  will  become  apparent  as  we  go 

along . 

The  trial  functions   are   distinguished  from  the   exact   functions  by  the 

subscript   t,    and  are   required  to  be   even  and  odd  and  to  have  the   asymptotic 

forms  given  by  equations    (9)   except   for  the   substitution   of  trial  phase 

shifts   n   ,    and  n  ^   for  the   exact  phase   shifts   -n     and  ri    .    In  addition,   the 
'et  'ot  'e  'o 

trial  function  is   required  to  be  normal  to  all   conducting  surfaces   and  both 
the   function  and  first   derivatives   are   required  to  be   continuous   in  the 
interior  of  the  waveguide. 
Consider  the   integral 


J  (E.X   E,   -    E^    -oCm^   --/(£•    VX  VX  E^   -   E,   • 


X  V  X  E)dT        (10) 


11  - 


where  the   range   of  Integration  is   over  the   interior  of  the  waveguide.      (it 
is   to  be  understood,    from  here  on^    that    'the   interior  of  the  waveguides' 
excludes   any  metallic  obstacles.)      The   second  tenn  in  the   integrand  on  the 
left  hand  side   of  equation   (lO)   vanishes  by  virtue   of  equation   (8a),    and 
the  right  hand  side   can  be  transformed  by  Green's   theorem  to  give 


j    E    •  (tC   E  dT  =     j    (5.  X  Y  X 


E      -   EX  V  X  E)    ?    dcy     . 


(11) 


The  surface  integral  vanishes  on  the  conducting  surfaces  since  E  and  E 
are  required  to  be  normal  to  these  surfaces,  and  the  contribution  from  the 
asymptotic  region  is 


j  Eg  •  ^Eg^dT  =  2k  cot(ri^0)  -  cot(ri-9) 


-.2 


e(x,y)  X  a 


do- 


(12) 


where  a  is  the  unit  vector  in  the  z  direction  and  the  surface  integral  is 
~z 

over  the  cross  section  of  the  wave  guide.   It  is  convenient  at  this  point 
to  adopt  a  normalization  and  notation  conforming  more  closely  to  that  of 
Kato.   Let  u  and  w  be  defined  by 


(13a) 


^9  =  ^0t  -  ^c 


(13b) 
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With  these  definitions,  equation  (l2)  can  be  vritten  as 

kcot(Ti-O)  =  kcot(Ti^-9)  -ij  Ug^-  dC  J^Q^^-^  -^  i   Ji£9*'^WQdT         (Ik) 

in  analogy  vith  equation  (3)«   The  factor  -^   arises  because  the  range  of 
integration  has  been  doubled.  The  first  tvo  tenns  on  the  right  hand  side 
of  equation (l4)  can  be  calculated  by  specifying  the  trial  function,  and 
constitute  a  variationaJ.  approximation  for  kcot(Ti-O).   The  third  term, 
z      ^Q  '         ^a^'^   '    ^^   *^^  error  term  and  is  of  the  order  of  the  square 
of  the  error  in  the  trial  function. 

We  now  derive  one  consequence  of  equation  (l^).   We  have  t^  as  the 
phase  shift  associated  with  W.   We  let  t]  +  d-r]  be  the  phase  shift  associated 
with  W -)-  dW.  where  dW(r)  >  0  for  all  r.   Applying  equation  (1^+)  to  the 
determination  of  tj  ^  dr],   we  choose  for  our  trial  function  the  exact  solu- 
tion u^  associated  with  W.   We  find 
—0 


dT^  =  (x  sin^Ti/k)   \^(^)    dW(r)  dr  >  0  , 


where  the  inequality  is  a  consequence  of  the  non-negativeness  of  dW(r)  and 
the  reality  of  u  (r).   More  generally,  it  then  follows  that  if  W  (r)  > 
W  (r),  then  'Hp  ^  It  •   This  result  is  generally  referred  to  as  the  mono- 
tonicity  theorem. 
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l^. .   Rigorous  bounds  on  the  error 

Kato's  method  for  obtaining  rigorous  bounds  on  the  error  term,  and  hence 
on  cot(T]-0).  can  be  applied  to  the  waveguide  case  with  minor  modification. 
In  order  to  accomplish  this,  we  make  use  of  the  following  information  about  w 


0* 


cCm^      =  ^<^  (15a) 


w^ 


C   e_(x,y)  cos(kz+9)   ,   z  -^  +  00  (15^) 


w^  ^»  C   e(x,y)  sin(kz+9)   ,   z  ->  +  00  (l5c) 

~0o       o   —  s    ^    I 


where  C   and  C   are  unknown  constants.  Also,  w_,  and  w^  are  even  and  odd 
e       o  — 9e     ~-9o 

functions  of  z  respectively;  w  is  noj:Tnal  to  all  conducting  surfaces,  and 

w  and  its  first  derivatives  are  continuous  in  the  interior  of  the  waveguide. 

Consider  the  equation 

c^0(r)  +  kip(r)^(r)   =0  (16) 

The  function  0(r)  is  required  to  be  normal  to  all  conducting  surfaces.   The 
weight  function  p(r)  >  0  is  an  even  function  of  z  ajid  tends  to  zero  at  large 
distances  from  the  obstacle  in  such  a  way  that    p(r)dT  converges.   Equation 
(16)  may  be  regarded  as  a  scattering  problem  of  the  type  of  equation  (8),  with  a 
'potential'  W  +  |i.p  .  The  even  and  odd  standing  wave  solutions  are  associated 
with  entirely  separate  problems  and  the  distinction  between  them  is  not  rele- 
vant to  the  following  discussion.   Let  5(|a)  be  the  phase  shift  corresponding 
to  this  'potential' .   (The  choice  of  p  must  be  such  that  this  problem  also 


Ik 


satisfies  the  four  restrictions  listed  in  the  introduction  in  order  that 

there  be  just  one  phase  shift  o(|a)  associated  with  each  value  of  [x.) 

Equation  \L 6} will  be  referred  to  as  the  associated  eigenvalue  problem  with 

eigenvalues  u  and  elsen   functions  0  (r)  when  ofu.  )  is  restricted  by  the 
t>         ^n  n  ~»         n  '' 

condition 

5(|a  )   =  9  +  n:rt   .  (l?) 

If  f_   and  g  are  two  functions  which  satisfy  the  boundary  conditions  appro- 
priate to  0  ,  then  it  follows  from  Green's  theorem  that 


I   f   •  oC   ^  ^^   =     j  ^  '  ^  t 


dT  (18) 


where  the  integrals  are  over  the  interior  of  the  waveguide.   Thus  the 
operator  c>C    restricted  by  these  boundary  conditions  is  hermitean;  it  follows 
that  the  eigenvalues  are  real  and  the  eigenfunctions  satisfy  orthogonality 
relations 


k   0  •  0  p  dT  =  5  (19) 

'  — m  ~-n  mn  ^ 


where  5   =  1  if  m  =  n  and  5   =  0  if  m  ^  n  and  the  normalization  is  chosen 

mn  mn 

for  convenience. 

Using  the  monotonicity  theorem  discussed  earlier,  it  follows  since 
p(r)  >  0  for  all  r  that 


^1^  >  0  .  (20) 

d|i    —  ^ 
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The  eigenvalues  therefore  satisfy  the  relationship 

^l  ^n   >   1^    .  (21) 

n+1  -   n  ^      ' 

It  can  be  shown  that  the  eigenfunctions  _^  form  a  complete  set  in  the 
sense  that  the  Parseval  identities  hold  for  any  functions  _F,  G  vhich  satisfy 
certain  symmetry  conditions  to  be  discussed  in  the  folloving  section,  which 
are  normal  to  all  conducting  surfaces  and  which  are  continuous  and  have  contin- 
uous first  derivatives,  that  is 

k  J  F^p  dT  =  _^  a^^  ,      k  J  F  .  Gp  dx  =  "^  a^b^  , 


k  j  G^p  dT  =  Yl     \' 


(22) 


the  a  and  b  are  Fourier  coefficients  defined  by 
n      n 


a  =  k 
n 


0^  .  Fp  dT  ,   b^  =  k  I  _0^  .  Gp  dT   .  (25) 


-1 


Wow  let  F  satisfy  the  boundary  conditions  (17)  and  set  G_  =  p   00  F. 
Then 

b   =kf0  -o^CFdT  =ke-^0  'Fdx  =-tikf0-FpdT 
n  J  ~n  —        J         J-n      —  "■  j  -- ^  — 

=  -   ^^n  (2h) 

where  we  have  made   use   of  the   hermltean  property  of  oC,    equation   (18). 
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The  Parseval  identities  become 


k  I  F^  pdT        =  X^  n^"^  b^^  (25a) 


k   F  •  o^ FdT      =-Yl     ki  '"^  ^  ^  (25^) 


k 


UFf  p-hx  =  i:  b  2  .  (25c) 


y 


n 


Let  a_  be  the  smallest  positive  eigenvalue  and  -P  the  smallest  (in  abso- 
9  '=' 

lute  value)  negative  eigenvalue,  so  that  -a    <  -[i^       <   Pq   •   Then  it 
follows  from  equations  (25)  that 

-a^"^  f  (^F)^  p"^dT  =  J  F  •  o<:  FdT  ^  Pq-^  J  {cClf   p"^dT   .      (26) 

The  point  is  that  while  c>C  is  an  unbounded  operator.  oC^        is  a  bounded 
operator,  its  eigenvalues  ranging  from  -P    to  a„   .   If  we  identify 
_F  with  w   defined  by  equation  (ijb)  and  recognize  that  06  w^  =  °^iig-(-^ 
equation  (26)  implies  that 

-^a^'^  J  Uu^^fp-hr  <    i[wQ'  ■/'w^dT  <^Pq-^  J  iaC  u^fp-hx    .   (2?) 

I  f  2   -1 

The  quantity    X  I  (0^-9+)   P~  '^'''  ^^  determined  by  specifying  the  trial 

function  J^^.-      Thus  the  problem  of  obtaining  rigorous  upper  and  lower  bounds 
on  the  error  term  in  equation  (l^)  is  reduced  to  that  of  obtaining  rigorous 
lower  bounds  on  a„  and  p   .  If  the  trial  function  u   is  a  good  approximation 
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to  the  exact  function  _u,  ^       ~o+^      ^^   ^  small  term  of  second  order,  and 
even  crude  lower  bounds  on  a  and  (3  can  provide  close  bounds  on  the  error 
.term  and  hence  on  cot(r|-9). 

5.    Completness 

The  question  of  completeness  Is  rather  more  complicated  for  the  present 
waveguide  case  than  for  the  one  independent  variable  quantum  mechanical  case 
discussed  by  Kato.   In  particular,  it  will  be  shown  that  the  0  as  defined 
by  the  differential  equation,  (16) .  and  by  the  boundary  condition, 
equation  (l7).  are  not  necessarily  complete.   Limitations  are  thereby  imposed 
on  the  trial  fimction,  in  addition  to  the  boundary  conditions  stated  in 
Section  5?  if  the  bounds  deduced  in  Section  h-   on  cot(T|-0)  are  to  be  valid. 
It  will  be  seen,  however,  that  the  limitations  are  in  no  sense  severe;  on 
the  contrary,  it  will  only  be  required  that  the  trial  function  have  the  same 
symmetry  properties  as  the  true  function.   Clearly,  no  restrictions  whatever 
are  thereby  introduced  on  the  accuracy  of  the  bounds  that  can  be  achieved. 
The  fact  remains  that  each  of  the  symmetry  properties  must  be  recognized  if 
the  Kato  bounds  are  to  be  valid. 

That  the  0  are  not  always  complete  can  be  seen  by  considering  the  simple 
example  of  a  waveguide  for  which  the  obstacle  Is  a  dielectric  slab  which 
extends  to  the  conducting  boundaries  of  the  waveguide  and  whose  permittivity 
is  Independent  of  x  and  y.   Then  each  of  the  0  will  everywhere  (not  only 
asymptotically)  have  the  same  x  and  y  dependence.   In  particular,  each  0  can 
be  expressed  as  a  product  of  e_(x,y)  and  a  function  of  z,  where  e(x,y)  is  the 
form  function  appropriate  to  the  dominant  mode,  and  the  problem  is  then 
effectively  one  dimensional.   In  this  example  the  0  are  obviously  not  a 
complete  set;   this  need  cause  no  concern,  however.   It  is  clear  that  the  true 
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solution  _u  would  also  have  the  same  x,y  dependence  as  the  dominant  mode. 
One  would  then  choose  the  trial  function  u  to  have  this  x,y  dependence,,  in 
which  case  w  =  u  -  u  would  also  have  this  x,y  dependence,  and  this  w 
could  be  expanded  in  terms  of  the  ^    .      On  the  other  hand,  if  one  failed  to 
recognize  the  one  dimensionality  of  the  problem  and  used  any  x,y  dependence 
other  than  that  given  by  e(x,y),,  w  could  not  be  expanded  in  terms  of  the  0 
and  the  bounds  obtained  on  cot(-ri-O)  in  Section  h   would  not  be  rigorous. 

In  discussing  the  question  of  completeness  for  an  arbitrary  W(x,y,z),  it 
is  assumed  for  simplicity  of  presentation  that  W  and  p  both  vanish  iden- 
tically beyond  |z|  =  d.   It  is  believed,  however,  that  the  conclusions  of 
this  section  can  be  established  under  more  general  circumstances.   Wow  it  can 
be  shown  that  the  eigenfunctions  of  a  positive  definite,  hermitian  operator 
form  a  complete  set.    Actually,  the  essential  point  is  not  that  the  operator 
be  positive  definite  but  rather  that  there  be  a  minimum  negative  eigenvalue. 
That  there  is  a  minimum  negative  eigenvalue  in  the  associated  eigenvalue 
problem  follows  from  the  assumption  that  W  and  p  vanish  identically  beyond 
|z|  =  d,  since  there  is  a  minimum  number  of  nodes  which  the  wave  functions 
can  have  in  the  region  -d  <  z  <  d. 

It  remains  to  determine  the  boundary  conditions  for  which  oC  is  hermitian. 
It  follows  from  Green's  theorem  that  for  two  functions  f(r)  and  g(r)  which 
are  continuous  and  have  continuous  first  derivatives,  and  which  are  normal 
to  the  conducting  boundaries  of  the  waveguide, 


•  M£  '  '"^  g  -  g  •  '=^£)dT  =  J(£XYXg-gxvxf) 


do- 


where  the  surface  integral  is  over  end  surfaces  of  the  waveguide  in  the 
asymptotic  region.   Consider  a  set  of  eigenfunctions  \|f  of  the  differential 
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equation 


— m    m   — m 


defined  by  the  boundary  condition 


1 


(i|f     X     VXiJf      -111      XVX\|r      )'dcr=0 
^  J.m  —        -n       -n         —         -m  — 


for  all  m  and  n.   Then  the  operator  o^    is  hermitiaji  with  respect  to  the  set  of 

functions  i|/  and  it  follovs  that  these  functions  form  a  complete  set.   The  set 

-m 

of  functions  -^     can  be  divided  into  tvo  subsets.   The  first  subset  consists 

— m 

of  functions  which  have  the  asymptotic  form  of  the  dominant  mode  and  whose 
phase  shifts  S(X  )  differ  by  integral  multiples  of  «•    This  subset  is  iden- 
tical with  a  set  of  associated  eigenfunctions  0  with  X  =  \J-   •      The  second 

-~xi  m    n 

subset  consists  of  functions  which  vanish  asymptotically.   These  functions  do 
not  occur  in  the  one  independent  variable  quantum  mechanical  problem  with 
fixed,  positive  energy,  and  it  follows  that  for  that  problem  each  of  the  t 
equals  one  of  the  jz5  which  then  form  a  complete  set.   However,  in  the  waveguide 
problem  the  functions  which  vanish  asymptotically  can  exist  under  special 
circumstances  which  will  now  be  considered. 

In  the  region  of  large  |z|,  where  W  and  p  are  zero,  each  4^  Is  a  super- 
position of  evanescent  waveguide  modes  and  of  the  dominant  waveguide  mode. 
(The  term  waveguide  mode  will  be  used  for  modes  in  the  absence  of  an  obstacle.) 
The  coefficient  of  the  dominant  waveguide  mode  will  be  denoted  by  C  .   There 
are  two  alternative  possibilities. 

fl)   If  each  C  is  different  from  zero,  each  t_  equals  one  of  the  0 
^  '  j]2  ■*-in.  —XI 

and  each  X  equals  one  of  the  \i   .      The  0  then  form  a  complete  set. 
m  n      ~n 
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(2)   If  C  =0  for  any  m,  the  corresponding  t  vanishes  asyinptotically 
and  the  0  do  not  then  form  a  complete  set. 

The  condition  C  =0  can  only  arise  if  some  of  the  evanescent  waveeuide 

m 

modes  are  not  coupled  to  the  dominant  waveguide  mode.   This  lack  of  coupling 
can  only  occur  if  W  and  p  have  some  special  symmetry  property.   We  have 
already  given  one  example  in  which  W  and  p  are  independent  of  the  transverse 
coordinates.   As  another  example,  we  consider  a  rectangular  waveguide  where 
the  z  axis  is  the  axis  of  the  waveguide,  the  y  axis  is  the  axis  of  polariza- 
tion of  the  dominant  mode  and  the  waveguide  extends  from  x  =  0  to  x  =  a. 
We  further  specify  that  W  and  p  are  independent  of  y  and  are  even  functions 

of  X  -  -ya.   Then  the  fimctions  0  share  the  polarization  of  the  dominant 
t  ~n 


mode,  are  independent  of  y,  and  are  even  functions  of  x  -  ^  a;  the  0^  are  not, 


therefore,  a  complete  set.   The  set  of  functions  ^     which  is  complete 

includes  functions  which  are  odd  functions  of  x  -  -£  a.,    functions  which 

depend  upon  y  and  whose  polarization  is  not  along  the  y  axis,  and  functions 

■\-rLth  both  of  these  properties,  all  of  which  vanish  asymptotically.   However, 

the  exact  solution  u  is  also  polarized  parallel  to  the  y  axis,  is  independent 

of  y,  and  is  an  even  function  of  x  -  g-  a.   If  the  trial  function  u  is  chosen 

to  have  the  same  properties  it  will  be  possible  to  expand  w  in  terms  of  the 

incomplete  set  0  .  .  • 

In  summary,  if  W  and  p  have  no  symmetry  properties,  the  "^     are  identical 

with  the  0  so  the  latter  form  a  complete  set.   On  the  other  hand,  if  W  and  p 
-n 

do  have  some  symmetry  properties,  the  functions  0  have  some  common  symmetry 
properties  and  do  not  then  form  a  complete  set.   The  set  of  functions  ^     is 
complete,  however,  since  it  then  contains  functions  which  vanish  asymptotic- 
ally and  have  different  symmetry  properties.   For  an  arbitrary  trial  function 
u,  the  difference  function  w  could  only  be  expanded  in  terms  of  the  J^^- 
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one  could  not  then  use  the  monotonicity  theorem  to  order  the  eigenvalues 
in  terms  of  the  phase  shifts,  an  ordering  which  is  essential  for  the  appli- 
cability of  the  method  used  in  the  following  section  for  the  determination 
of  bounds  on  a„  and  p  .   It  is  of  course  quite  possible  that  methods  could  be 
devised  for  the  determination  of  bounds  on  a_  and  p  even  for  the  case  for 
which  evanescent  modes  must  be  included;  if,  however,  ii  is  chosen  with  the 
correct  symmetry,  w  can  be  expanded  in  terms  of  the  0  and  the  question  does 
not  then  arise.  ■  v - 


6.   Lower  bounds  on  a^  and  p  ■■  ''.  '  ■    '-^  ?2,..':'  ) 

The  general  procedure  for  obtaining  lower  bounds  on  a^  ajid  p  involves 
the  use  of  comparison  potentials  for  which  the  scattering  problem  can  be 
solved  exactly,  and  the  monotonicity  theorem,  equation  (20).  A  nvraiber  of 
cases  have  been  treated  in  the  literature.  '    In  this  section,  we  will 
consider  just  three  special  cases  applicable  to  a  wide  variety  of  waveguide 
scattering  problems.  We  restrict  ourselves  to  waveguides  which  are  of  uni- 
form cross-section  everywhere.   These  three  cases  are  particularly  simple 
to  apply. 

In  the  context  of  the  following  discussion,  metallic  obstacles  may  be 
regarded  as  dielectric  obstacles  with  infinite  negative  permittivity.   (This 
comment  does  not  apply  to  the  preceding  discussion  of  the  variational  prin- 
ciple and  associated  eigenvalue  because  the  integrals  would  be 
divergent  with  an  infinite  'potential'.)   Suppose  the  obstacle  does  not 
extend  beyond  z  =  +  d,  and  that  p  is  identically  zero  for  [zj  >  d  and 
positive  for  |z|  <  d. 

We  will  first  obtain  a  lower  bound  on  p^.   As  p.  approaches  -  oo,  it  is 
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evident  that  the  phase  shift  5(|_i)  of  the  associated  eigenvalue  problem 
approaches  -kd  in  the  odd  case  and  -kd  -  3  jt  in  the  even  case.   Now  if 
-kd  >  Q  -  jt  in  the  odd  case  or  -  kd  -  j  jt  >  9  -  jt  in  the  even  case  and  if 
T|  <  0,  then  as  shown  in  Figure  1;,  there  can  be  no  negative  eigenvalues  and 
P  can  be  regarded  as  infinite.   In  this  case  the  error  term  in  equation 
(1^)  is  always  negative,  and  the  variational  approximation  to  k  cot(T]-0) 
is  always  an  upper  bound  on  the  exact  value. 

A  lower  bound  on  a_  can  be  obtained  as  follows.  We  know  from  the  mono- 
tonicity  theorem,  equation  (20),  that  the  phase  shift  6(fi)  is  less  than  the 
phase  shift  6'(|j.)  which  would  result  if  the  obstacle  were  replaced  by  a 
dielectric  slab  which  fills  the  space  -d  <  z  <  d  out  to  the  conducting 
boundaries  of  the  waveguide,  and  whose  'potential'  is  equal  to  the  maximum 
value  of  ¥  +  p.p.   This  phase  shift  is  given  by 


5'  (n)  =  -kd  +  tan"-^(  J^  tan  k'd)  (28a) 

e  ^  kd 


5'  (li)  =  -kd  +  tan  ^[   p^  tan  k'd)  (28b) 


k'  =  k'(u)  =  [k^  +  (W  +  Mp)   1^  (28c) 


It  can  be  seen  from  Figure  2  that  a  useful  lower  bound  on  a^  is  provided 
by  the  value  of  n,  to  be  denoted  by  a^',  defined  by  5' (a')  =   0,  provided  that 
T)'  =  5'(o)  <  0  and  ri  >  0  -rt  .   For  a  dielectric  obstacle  t]  will  always  be 
greater  than  zero  and  the  last  condition  will  automatically  be  satisfied. 
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e-TT 


Figure   1 


Determination  of  a  lower  bound  on  p.  assuming  Vl  <  0,  and 

assuming  -kd  >  0  -  jt  in  the  odd  case  and  -kd  -  -£«  >  0  -  jt 

in  the  even  case.   Since  6(|-i)  must  exceed  0  -  tc,  P   is  infinite. 
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~e-ir 


Figure   2 


Determination  of  a  lower  bound,  on  a  assuming  "1^ '  <  0  and 
1^  >  0  -  jt.   The  upper  curve  corresponds  to  a  dielectric  slab 
filling  the  guide  in  the  range  -d  <  z  <  d,  whose  'potential' 
equals  the  maximum  'potential'  of  the  obstacle.   From  the  mono- 
tonicity  theorem,  6'(|j)  >  5(ji)  and  therefore  a  '  <  a   . 
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From  the  rough  bounds  already  obtained  on  t], 


-kd  --Jjt  <  Tig  <  T]^  (29a) 


-kd  <  n  <  ri'  (29b) 

'o    'o 


it  follows  that  a  sufficient  condition  for  obtaining  Ic-.'sr  bounds  on  a^  and 

2        -^ 
6„  by  the  procedure  described  above  is  that  d(k  +  W   )   be  less  than  4  rt 
^0  "^     -^  max  * 

in  the  even  case  and  less  than  «  in  the  odd  case.  An  equivalent  statement 
is  that  2d  must  be  less  than  y  X  in  the  even  case  and  X  in  the  odd  case, 
vhere  X  is  the  smallest  guide  wavelength  in  the  dielectric.   However,  it 
should  be  noted  that  these  restrictions  are  by  no  means  fundamental,  since 
there  are  many  more  possibilities  for  obtaining  lower  bounds  on  a_  and  Pg. 

There  is  one  case  of  interest  in  which  p   can  be  determined  exactly. 
If  W  is  everywhere  positive  we  can  choose  p  =  W.   Then  the  'potential'  of 
the  associated  eigenvalue,  W  +  \ip.    vanishes  when  n  =  -  1  and  the  phase  shift 
5(-l)  vanishes  also.   If  we  choose  0=0,  it  is  evident  that  p.  =  -1  is  an 
eigenvalue.   Furthermore,  if  we  know  that  x]  <  n ,    it  is  the  smallest 
(in  absolute  value)  negative  eigenvalue  as  can  be  seen  from  Figure  3.  Thus 
for  the  case  in  which  W>0,  p=W,  0=0  and  ri  <  n,  the  eigenvalue  p 
is  given  exactly  by 


P   =  1  (30) 


o 


(More  generally,  if  W  >  0,  if  p  =  W,  and  if  ti  <  0,  then  p^  >  1.) 
Equation  (30)  is  applied  in  Part  2  to  a  problem  in  which  a  special  trial 
function  is  used  which  generates  the  Schwinger  integral  variational  principle. 
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Figure  3 


Determination  of  the  exact  value  of  p  assuming  W  >  0,  p  =  W^ 
9  =  0  and  T\^  <  «.   In  this  case,  p   is  given  exactly  by  p   =  1. 
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Lambert  Saint-Louis  Municipal  Airport 

Box  516,   St.  Louis  3,   Missouri 

Attn:     R.  D.  Detrlch,   Engineering  Library 

McMillan  Laboratory,    Inc. 

Prownvllle  Avenue 

Ipswich,   Massachusetts 

Attn:     Security  Officer,    Document  Room 

Melpar,  Inc.  (2) 

30O0  Arlington  Boulevard 

Falls  Church,  Virginia 

Attn:  Engineering  Technical  Library 

Microwave  Development  Laboratory 
90  Broad  Street 
Babson  Park  57,  Massachusetts 
Attn:  N.  Tucker,  General  Manager 

Microwave  Radiation  Company  Inc. 

19223  South  Hamilton  Street 

Oardena,  California 

Attn:  Mr.  Morris  J.  Ehrlich,  President 

Chance  Vought  Aircraft,  Ino. 
9311j  West  Jefferson  Street 
Dallas,  Texas 
Attn:  Mr.  H.  S.  White,  Librarian 

Northrop  Aircraft,  Inc. 
Hawthorne,  California 

Attn:  Mr.  E.  A.  Freltas,  Library  Deptjll:5 
inoi  E.  Broadway 

Remington  Rand  Unly,  -  Division  of  Sperry 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,  Pennsylvania 
Attn:  Mr.  John  F.  McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,  Inc. 
12211:  Lakewood  Boulevard 
Downey,  California 
Attn:  Engineering  Library  1)95-115 

North  American  Aviation,  Inc. 
Los  Angeles  International  Airport 
Los  Angeles  1:5,  California 
Attn:  Engineering  Technical  File 

Page  Communications  Engineers,  Inc. 
710  Fourteenth  Street,  Northwest 
Washington  5,  D.  C. 
Attn:  Librarian 

Philco  Corporation  Research  Division 

Branch  Library 

ti700  Wissachickon  Avenue 

Philadelphia  hk.   Pa. 

Attn:  Mrs.  Dorothy  S.  Collins 


Pickard  and  Bums,  Inc. 

2llO  Highland  Avenue 

Needhan  9b,  Mass. 

Attn:  Dr.  J.  T.  DeBettenoourt 

Polytechnic  Research  and  Development 

Company,  Inc. 

202  Tlllary  Street 

Brooklyn  1,  New  Tork 

Attn:  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 
Attn:  Dr.  John  Ruze 

Radiation,  Inc. 
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Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.L.  Cox 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  Tork 

Attn:  P.  S.  Carter,  Lab.  Library 

RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attn:  Miss  Fern  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Administrator 

The  Ramo-Wooldridge  Corporation 
P.O.  Box  1:51)53  Airport  Station 
Los  Angeles  Ii5,  California 
Attn:  Margaret  C.  Whitnah, 
Chief  Librarian 

Hoover  Microwave  Co. 
9592  Baltimore  Avenue 
College  Park,  Maryland 

Director,  HSAF  Project  RAND 

7iai  Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:  Grace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Raytheon  Manufacturing  Companv 
Wayland  Laboratory,  State  Road 
Wayland,  Mass . 
Attn:  Mr.  Robert  Borts 

Raytheon  Manufacturing  Comrany 
Wayland  Laboratory 
Wayland,  Mass. 

Attn:  Miss  Alice  G,  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmlngdale,  Long  Island,  N,  Y. 
Attn:  Engineering  Library 

Thru:  Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmlngdale,  Long  Island,  N,Y. 

Rheen  Manufacturing  Company 
9236  East  Hall  Road 
Downey,  California 
Attn:  J.  C.  Joerger 

Trans-Tech,  Inc. 
P.  0.  i?ox  31:6 
Frederick,  Maryland 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  California 
Attn:  Librae/  -  unclassified 

Sage  Laboratories 
159  Linden  Street 
Wellesley  81,  Mass. 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Library 

Sandia  Corporation,  Sandla  Base 

P.O.  Box  5900,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  W.  Tumbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attn:  Library,  Engineering  Division 

Sylvanla  Electric  Products,  Inc. 
100  First  Avenue 
Waltham  51),  Mass. 

Attn:  Charles  A.  Thornhill,  Report  Librarian 
Walthan  Laboratories  Library 

Systems  Laboratories  Corporation 
li)''52  Ventura  Boulevard 
Sherman  Oaks,  California 
Attn:  Donald  L.  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:  M.  L.  Henderson,  Librarian 

A.  S.  Thomas,  Inc. 

I6l  Devonshire  Street 

Boston  10,  Mass. 

Attn:  A.  3.  Thomas,  President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 

Chu  Associates 
P.  0.  Box  3fi7 
Whitcomb  Avenue 
Littleton,  Mass. 

Microwave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Companv 
Missile  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  501),  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Command 

l!)th  Street  and  Constitution  Avenue 

Washington,   D.  C, 

Attn:     Mr.  Ottlng,   SRY 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l  Airport  Box  71)6 
Baltimore  3,   Maryland 
Attn:     Engineering  Library 
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Wheeler  Laboratories,  Inc, 
122  Catter  Mill  Road 
Oreat  Neck,  New  York 
Attn:  Mr.  Harold  A,  Wheeler 

Zenith  Plastics  Go. 
Box  91 

Gardena,  California 
Attn:  Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute 
of  the  University  of  Alaska 
College 
Alaska 

Dniverslty  of  California 

Berkeley  U,  California 

Attn:  Dr.  Samuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec.  Eng.  Electronics 
Research  Lab, 

University  of  California 
Electronics  Research  Lab. 
332  Cory  Hall 
Berkeley  It,  California 
Attn:  J.  R.  Whinnery 

Callfomia  Institute  of  Technology 
Jet  Propulsion  Laboratory 
IjPOO  Oak  Grove  Drive 
Pasadena,  California 
Attn:  Mr.  I.  E.  Mewlan 

California  Institute  of  Technology 
1201  E.  Calif  omi  a  Street 
Pasadena,  California 
Attn:  Dr.  C.  Papas 

Carnegie  Institute  of  Technology, 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Prof.  A.  E.  Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.  C.  Dalnan 

University  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,  Florida 

Attn:     Prof.  M.   H.  Latour,   Library 

Library 

Georgia  Institute  of  Technology 

Engineering  Experiment  Station 

Atlanta,  Georgia 

Attn:  Mrs.  J.H.  Crosland,  Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,  303A  Pierce  Hall 

Oxford  Street,  Cambridge  3",  Mass. 

Attn:  Mrs.  E.L.  Hufschmidt,  Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,  Mass. 
Attn:  Dr.  Fred  L.  Whipple 

University  of  Illinois 
Documents  Division  Library 
Urbana,  Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,  Illinois 

Attn:  Dr.  P.  E.  Moyes,  Department  of 
Electrical  Engineering 

The  Johns  HopAcins  University 
Homewood  Campus 
Department  of  Physics 
Baltimore  18,  Maryland 
Attn:  Dr.  Donald  E.  Kerr 

Sandia  Corporation 
Attn:  Organization  1U23 
Sandia  Base 
Albuquerque,  New  Mexico 


Applied  Physics  Laboratory 
The  JoWPHopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Mr.  George  L.  Selelstad 

Massachusetts  Institute  of  Technology 

Research  Laboratory  of  Electronics 

Room  20B-221 

Cambridge  39,  Haasachuaetts 

Attn:  John  H.  Hewitt 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Mass. 

Attn:  Document  Room  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J.  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:  Prof.  K.  M.  Siegel 
912  N.  Main  St., 
Ann  Arbor,  Michigan 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Tpailanti,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  lit,  Minnesota 

Attn:  Mr.  Robert  H.  Stuiran,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  R.  E,  Bean 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,  Ohio 
Attn:  Dr.  T.E.  Tice 

Dept.  of  Elec.  Engineering 

The  UnlTersity  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attn:  Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Broolflyn 

Microwave  Research  Institute 

??  Johnson  Street 

Brookljm,  New  York 

Attn:  Dr.  Arthur  A.  Ollner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
5$  Johnson  Street 
Brooklyn,  New  York 
Attn:  Mr.  A.  E.  Laemmel 

Syracuse  University  Research  Institute 
Collendale  Canpua 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  0.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Gerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W,  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn:  Prof.  0.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R.  Mingins 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G.  Held,  Associate  Professor 

Stanford  University 
Stanford,  California 
Attn:  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  laboratory 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:  Mr.  H.  W.  Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Dr.  C.  M.  Angulo 

Case  Institute  of  Technology 

Cleveland,  Ohio 

Attn:  Prof.  S,  Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momingside  Heights 

New  York,  N.  Y. 

Attn:  Dr.  Schlesinger 

McGill  University 

Montreal,  Canada 

Attn:  Prof.  0.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schulta 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Peimsylvania 

University  of  Pennsylvania 

Institute  of  Cooperative  Research 

3U00  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W.  Cumberland  Avenue 
Knoxville  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Depai^ment  of  Electrical  Engineering 

Seattle,  Washington 

Attn:  Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Attn:  Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C.  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  h,    Pennsylvania 


Ionosphere  Research  Laboratory 

Pennsylvania  State  College 

State  College,  Pennsylvania 

fTTNl  Professor  A.  H.  Wayniolt,  Director 


Institute  of  Mathematical  Sciences 
25  Wairerly  Place 
New  Tork  3,  New  York 
ATTNt  Librarian 
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Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
ATTNi  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.  C. 

ATTN:  Dr.  W.  K.  Saunders 

Applied  Mathematics  and  Statistic-s  Lab. 

Stanford  University 

Stanford,  aiifomia 

ATTN;  Dr.  Albert  H.  Bowker 

Department  of  Physics  and  Astronomy 
Michigan  State  College 
East  Lansing,  Michigan 
ATTNi  Dr.  A.  Leitner 

Universitv  of  Tennessee 
Knox'rtlle,  Tennessee 
ATTN:  Dr.  Fred  A.  Fioken 

Lebanon  Valley  College 
Annville,  Pennsylvania 
ATTN:  Professor  B.H.  Bissinger 

General  Atomic 

P.  0.  Box  608 

San  Diego  12,  California 

Attn:  Dr.  Edward  Gerjuoy 

Department  of  Physics 
Amherst  College 
Amherst,  Mass. 
ATTN:  Dr.  Arnold  Arons 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
ATTN:  Dr.  A.  Erdelyi 

Mathematics  Department 

Stanford  University 
Stanford,  California 
ATTN;  Dr.  Harold  Levins 

University  of  Minnesota 

Minnpapolis  I'l,  Minnesota 

ATTN;  Professor  Paul  C.  Rosenbloom 

Department  of  Mathematics 

Stanford  University 

Stanford,  California 

ATTN:  Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  Johrf Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring.  Maryland 
ATTN:  Dr.  B.  S.  Gourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,  D.  C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Mass. 
ATTN:  Dr.  L.  J.  Chu 

Nuclear  Developnent  Associates,  Inc. 

5  New  Street 

White  Plains,  New  York 

ATTN;  Library 

California  Institute  of  Technology 
Electrical  Engineering 
Pasadena,  California 
ATTN;  Dr.  Zohrab  A.  Kaprielian 


Dr.  Rodman  Doll 
311  W.  Cross  Street 
Ypsilanti,  Michigan 

California  Inst,  of  Technology 
Pasadena,  California 
ATTN:  Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 
Woods  Hole  Oceanographio  Institute 
Woods  Hole,  Mass- 
National  Bureau  of  Standards 
Boulder,  Colorado 
ATTN:  Dr.  R.  Gallet 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  California 

University  of  Minnesota 
The  University  of  Library 
Minneapolis  lb,  Minnesota 
ATTN;  Exchange  Division 

Department  of  Mathematics 
University  of  California 
Berkeley,  California 
ATTN:  Professor  Bernard  Friedman 

Lincoln  Laboratory 

Massachusetts  Institute  of  Technology 

P.  0.  Box  73 

Lexington  73,  Massachusetts 

ATTN;  Dr.  Shou  Chin  Wang,  Room  C-351 

Melpar,  Inc., 

3000  Arlington  Boulevard 

falls  Church,  Virginia 

ATTN:  Mr.  K.  S.  Kelleher,  Section  Head 


Hq.  Air  Force  Cambridge  Research  Center 

Laurence  0.  Hanscom  Field 

Bedford,  Mass. 

ATTN;  Mr.  Francis  J.  Zucker,  CRRD 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTN;  Dr.  Philip  Newman,  CRRK 

Mr.  N.  0.  Gerson 

Trapelo  Road 

South  Lincoln,  Mass. 

Dr.  Richard  B.  Barrar 
Systems  Development  Corp. 
2li00  Colorado  Avenue 
Santa  Monica,  California 

Columbia  University  Hudson  Laboratories 
P.O.  Box  239 

1U5  Palisade  Street,  Dobbs  Ferry,  N.  Y. 
ATTN;  Dr.  N.  W.  Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
University  of  Maryland 
College  Park,  Maryland 
ATTN;  Dr.  Elliott  Montroll 

Department  of  Electrical  Engineering 

Washington  University 

Saint  Louis  5,  Mo. 

ATTN:  Professor  J.  Van  Bladel 

Department  of  the  Navy 

Office  of  Naval  Research  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,  California 

Brandeis  Tlniversity 
Waltham,  Mass. 
ATTN:  Library 

General  Electric  Company 
Mlctowave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,  California 
ATTN:  Library 


Smyth  Research  Associates 
3555  Aero  Court 
San  Diego  3,  California 
ATTN;  Dr.  John  B.  Smyth 

Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  California 
ATTN:  DrJjeorges  G.  Weill 

Naval  Research  Laboratory 

Washington  25,  D.  C. 

ATTN;  Henry  J.  Passerini,  Code  5Z75A 

Dr.  George  Kear 
5  Culver  Court 
Orinda,  California 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,  New  York 

ATTN;  Dr.  Nathan  Marcuvitz 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Jerry  Shmoys 

Department  of  Mathematics 
University  of  New  Mexico 
Albuquerque,  New  Mexico 
ATTN:  Dr.  I.  Kolodner 

Mathematics  Department 
Polytechnic  Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,  New  York 
ATTN;  Dr.  Harry  Hochstadt 

Ballistics  Research  Laboratory 
Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 
ATTN:  Dr.  PuUen  Keats 


Dr.  Lester  Kraus 
lt935  Whitehaven  Way 
San  Diego,  California 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis,  Minnesota 
Attn:  Dean  Athelston  Spllhaus 

Ohio  State  University 
Columbus,  Ohio 
Attn:  Prof.  C.  T.  Tai 
Department  of  Electrical  Eng. 

Naval  Research  Laboratories 

Washington  25,  D.  C. 

Attn:  W.  S.  Ament,  Code  5271 

Naval  Research  Laboratory 
Washington  25.  D.  C. 
Attn:  Dr.  Leslie  G.  McCracken,  Jr. 
Code  3933A 

Office  of  Naval  Research 
Deoartment  of  the  Navy 
Attn:  Geophysics  Branch,  Code  Ul6 
Washington  25,  D.  C, 

Office  of  Chief  Signal  Officer 
Signal  Plans  and  Operations  Division 
Attn:  SIOOL-2,  Room  20 
Cora.  Liaison  Br.,  Radio  Prop.  Sect. 
Th(?  Pentagon,  Washington  25,  D.  C. 

Defence  Research  Member 
Canadian  Joint  Staff 
2001  Connecticut  Street 
Washington,  D.  C. 

Central  Radio  Prop.  Lab. 
National  Bureau  of  Standards 
Attn:  Tochnioal  Reports  Library 
Boulder,  Colorado 

U.  S.  Weather  Bureau 
U.  S.  Department  of  Commerce 
Washington  25,  D.  C. 
Attn:  Dr.  Harry  Wexler 
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Federal  Communications  Commission 

Washin.rton  25,  D.  C. 

Attni  Mrs.  Barbara  C.  Grimes,  Librarian 

Upper  Atmosphere  Research  Section 
Central  Radio  Propagation  Laboratory 
National  Bureau  of  Standards 
Boulder,  Colorado 

Argonne  N'ational  Laboratory 

P.O.  Box  299 

Lemont ,    nilnols 

Attn:     Dr.   Hoylande  D.  Young 


Bell  Telephone  Labs. 
Murray  Hill,   New  Jersey 
Attn:     Dr.  S.   0.   Rice,   3B 
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Carnegie  Institute  of  Washington 
Dept.  of  Terrestrial  Magnetism 
S2lil  Broad  Branch  Road,  N.  W. 
Washington  15,  D.  C. 
Attn:  Library 

Georgia  Tech  Research  Institute 
225  N.  Avenue,  N.  W. 
Attn:  Dr.  James  E.  Boyd 
Atlanta,  Georgia 

University  of  Maryland 
College  Park,  Maryland 
Attn;  Dr.  A.  Weinstein 
Institute  of  Fluid  Dynamics 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Massachusetts 

Attn:  Prof.  Radford,  Division  3  Head 

Willow  Run  Research  Center 
University  of  Michigan 
Willow  Run  Airport 
Tpsilanti,  Michigan 
Attn;  Dr.  C.  L.  Dolph 

School  of  Engineering 
New  York  University 
University  Heights 
New  York,  New  York 

Shell  Fellowship  Committee  of  the 
Shell  Companies  Foundation,  Inc. 
SO  West  50th  Street 
New  York  20,  N.  Y. 
Attn:  Mr.  J.  R.  Janssen 

Esso  Research  and  Engineering  Co, 

P.  0.  Box  51 

Linden,    New  Jersey 

Attn;     Mr.  C.  L.  Brown,   Manager 

Union  Carbide  and  Carbon  Corp. 
30  E.   Ii2nd  Street 
New  York  17,   New  York 
Attn:     Mr.  L.   E.   Erlandson 
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